‘TREX source 


The length of a Berier curve is an elliptic integral. This is a Bezier curve. 


x =agt? +agt? +a,t+ 29 


y = bgt? + byt? + bit + yo 


The control points of a Bezier curve are: 
(r0,Yo) (t1,y1) (2,2) (x3, ys) 
This is the curve for those control points. 


a, = 3(2#1 — Xo) a2 = 3 (a — 2214+ £2) a3 = —% +321 -—32%9+ 23 


b,=3(y1—yo) b2 =3(yo—2y1+y2) 03 =—yot3y1 —3y2+ y3 


The length of the Bezier cuvre form 0 to u is: 


| V (ay + 2agt + 3.43 t?)? + (bj + 2bot + 3 bg t?)? dt 
0 
The control points of an example Bezier curve are: 
(to =6,y0=8) (1=1,y1=10) (@2=7,y2=3) (@3=4, 43 = 4) 


Solve for ay, a2, 43, by, bo, b3. 


a, = —15 ag = 33 a3 = 20 b; =6 bg = 27 b3 = 17 
The Bezier curve for this example is: 
eg = —200° +3347 -15t+6 y=17t2—2727+6t4+8 


For this example the length is: 


1 
3 | V 689 t4 — 1492 t3 + 1076 t? — 292% + 29 dt 
0 


Numerical integration gives 7.237223368328592885826619210956022413967067986017504163362886516 
DERIVE factored the polynomial: 


689 t* — 1492 t? + 10764? — 202¢ + 29 = 


Vv 23377 Vv 2 7 
eso (44 OUP apis 6005 23377 Fingal 
1378 474721 689 1378 474721 689 

i V6005 “i 23377 373 V6005 23377 7 

1378 474721 689 1378 474721 689 

ea V6005 i 23377 373 ve V6005 23377 7 

1378 474721 689 1378 474721 689 


V6005 23377 
1378 474721 


V6005 2 23377 373 7 oe 
1378 474721 689 689 


Now define values k,m, a, b: 
k=3V689 m=(1,1,1,1]) b=[1,1,1,1] e = [e1,e2,e3, e4] 
a = [—0.216589 + 0.0937743 7, —0.216589 — 0.0937743 2, —0.866139 + 0.0734550 7, —0.866139 — 0.0734550 7] 


Do a partial fraction expansion. 
I(m) = D219(1,m, 4, a, b, e) 


I(m) = I(4e1) + (1.29909 + 0.375097 2) I(3 e) + (0.383342 + 0.365466 2) I(2e1) 
—(0.0228475 — 0.0784923 2) I(e1) 


Reduce [(4 1). 
R35(0, 1,4, a, b,e) = 0 


A(2e1 +e2 +e3 + e4) + 31 (41) + (8.24774 + 0.937743 2) (3 e1) + (0.766684 + 0.730933 2) I(2 e1) 
—(0.0342712 — 0.117387) I(e1) = 0 
A(2e1 + eg + e3 + e4) = —0.138815 + 0.00794130% 
Substitute [(4 e). 
I(m) = (0.216516 + 0.0625162 7) I(3 e1) + (0.127780 + 0.121822 7) I(2e1) 


—(0.0114237 — 0.0392461 7) I(e,) + 0.0462716 — 0.00264710% 


Reduce (3 1). 
R35(—1,1,4,a, b, e) = 0 


A(ex +eg+e3+ €4) +2 I(3 €1) + (1.94864 + 0.562645 i) I(2 €1) 
+ (0.383342 + 0.365466 2) I(e1) — (0.0114237 — 0.0392461 7) T(0) = 0 
A(ez + eg + e3 + e4) = —0.0846852 


Substitute [(3 e1). 
I(m) = —0.0655894 I(2 e,) — (0.0415000 + 0.0123012 7) Tex) 


+(0.00246348 — 0.00389163 7) T(0) + 0.0554395 


Reduce [(2 1). 
R35(—2,1,4, a, b,e) = 0 


A(eg + e3 + e4) + I(2e1) + (0.0114237 — 0.0392461 7) T(—e1) 
+(0.649549 + 0.187548 2) I(e1) =0 
A(eg + e3 + e4) = —0.949300 — 0.327220% 
Substitute [(2e;). 
I(m) = (0.000749279 — 0.00257413 7) I(—e,) + 0.00110358 I(e1) 

+(0.00246348 — 0.00389163 7) T(0) — 0.00682448 — 0.0214622% 

Calculate the basic integrals numerically. 
I(—e,) = —17.2709 + 33.51427 I(0) = 12.0414 T(e,) = —3.86310 — 1.129177 
Substitute the values of the integrals. 
I(m) = 0.0919054 k I(m) = 7.23722 


That is a good match. 
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